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Abstract 

We consider the quantum state describing the Disoriented Chiral Condensate (DCC) 
which may be produced in high energy colhsions. Using the approach suggested by Ra- 
jagopal and Wilczek to describe the amphfication of the long wavelength classical pion 
modes, we consider the quantum-mechanical evolution of the initial vacuum state into the 
final squeezed state describing the DCC. The obtained wave function has some interesting 
properties which are discussed. 
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The possibility of the Disoriented Chiral Condensate (DCC) production in high energy 



hadronic or heavy ion colhsions attracted a lot of attention in the last years |]T|-||T4 . 

It is well known that QCD Lagrangian is invariant (approximately iif nonzero masses 
for the light Nj quarks are taken into account) under global chiral SU^Nf)^ x SU{Nf)ji, 
where Nj is the number of the light flavours. This symmetry is spontaneously broken down 
to vector SU{Nf)v which leads to Nj — 1 (quasi)goldstone bosons - pions (if Nf = 2) 
or pions, kaons and t] meson (if Nf = 3). The order parameter for this breaking is 
the vacuum expectation value of quark condensate < ipip >■ However one can imagine 
that under some special conditions, for example after high-energy collision, there is a 
"cool" region surrounded by a "hot" relatively thin expanding shell, which separates the 
internal region from the outer space. This picture was suggested by Bjorken and is called 
now the "Baked Alaska" scenario. In result the quark condensate may be disoriented in 
isotopical space - one gets the DCC. After hadronization the interior disoriented vacuum 
will collapse decaying into pions. The interesting signature of this process will be the 
coherent production of either charged or neutral pions. There are some arguments that 



the DCC was observed in the Centauro cosmic ray events |]T5 . 

It is convenient to consider the toy model describing the chiral dynamics - the linear 
sigma- model with four component field 0" = (cr, vf), where a and vr are an isoscalar and 
an isovector fields (here we use the same notation as in [^) 

s = I d^x[^d,rd^.r - ^irr - vY + Ha] (i) 

where H ~ describes the small explicit chiral symmetry breaking. The pion mass is 
ml = H/ f^, = \{fl — f^), where /^r =< cr >. The sigma mass is mo- = 2A/^. In strong 
coupling limit A ^ oo one gets the constraint + tt^ = In the usual vacuum one 
has <cr>=/^, <7f>=0 and since a is an isoscalar there is an isoscalar condensate 

< ipip > only. However one can consider another configuration - < a >= /7rCos6' and 

< 7? >= /jrnsin^, here n is some unit vector in the isospace, which describes DCC, i.e. 
some classical pion field configuration, which is metastable and decays after some time 
into pions - the signature for this event will be the large number of either neutral (vr^) or 
charged (vr^) pions. Using the classical picture of DCC which predicts equal probability 
for all isotopical orientation of condensate one can get [0] - @] the probability l/v7 for 
the fraction of neutral pion / = N.„o/{Npi+ + N^,- + A^.^o). 

The question we would like to discuss here is connected with the quantum description 
of the DCC . The simplest possibility is to describe it by the usual coherent state. However 
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the usual coherent state wave function was critisised in 0] because such a description 
leads to creation of pion system with arbitrary large charge fluctuations. Instead it was 
suggested in |^ that the quantum state must be isosinglet. For the state with 2n neutral 
pions and total number of pions 2N it was suggested: 

|^>~(2ataL- (4)2)^10 > (2) 

Then the probability to have 2n neutral pions is ||T6|, 0] 



^'"■^>= (2iV+l)!(y)^ ~V^- 

and corresponds to l/^ff distribution in classical picture. 

However let us note that the same distribution will be obtained and in the case of 
the arbitrary relative factor between charged and neutral creation operators [2a^|_al — 
exp(i^)(ao)^]^ and one can see that from this point of view the zero isospin condition is 
not so important - what was important indeed is that the wave function was constructed 
by the operators quadratic in and by the construction there are equal numbers of 7r+ 
and 7r~. 

It is interesting to understand what are the most natural class of these functions and 
what are the dynamical mechanisms leading to the generation of these functions. As 
we shall demonstrate here the quantum state of the DCC is the squeezed state. These 
quantum states known for a long time in quantum optics and measurement theory ( for 
a review of squeezed states see, for example -0)- The simplest one-mode squeezed 
state is parametrized by the two parameters r and and can be obtained by acting on 
the vacuum by the squeezed operators S'(r, 0) 

5(r,0)|O >= exp[^(e-''V - e^'^a^'')]\Q > (4) 

These states are minimum uncertainty states with AX1AX2 = 1/4, where a(a"l') = 
Xi ± X2, as well as coherent states exp(aa''")|0 >. However coherent states have some 
minimal quantum noise and AXi = AX2 = 1/2. For the squeezed states one can 
reduce the quantum noise for one variable (increasing it for the conjugate one) and 
AYi = e~''"/2, Ay2 = e''"/2 where Yi + iY2 = {Xi + 1X2)6'^'^. The mean particle number 
< >= sinh^ r. 

To get the squeezed quantum state for the DCC let us consider the mechanism of the 
amplification of the long wavelength pion modes, suggested by Rajagopal and Wilczek in 
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paper [P, where the the dynamics of the 0(4) hnear sigma model after quenching was 
considered. The amphfication of the long wavelength pion modes was found in the period 
immediately after quenching. This amplification leads to the coherent pion oscillations, 
i.e. to the creation of the DCC . Such a behaviour can be understood if one considers the 
equation of motion for the pion field OB 



d_ 

dt 



^7r{k, t) + [e + A(< <P' > (t) - v')]n{k, t) = (5) 



where we substituted the in the nonlinear term in (|]) by its spatial average < 0^ > 
(t) - this is nothing but the Hartree-Fock or mean field approximation. In the initial 
conditions one has < cfP' >< v"^ and the long wavelength modes of the pion field with 
fc^ < A(f^— < 0^ >) start growing exponentially. The < cfP' > (t) starts to oscillate near 
the vacuum expectation value < a > and after some time the oscillations will be damped 
enough so that all the modes will be stable. Thus we see that at classical level each long 
wavelength mode is described by the equation for a parametrically excited oscillator and 
one get the DCC in result of amplification of the zero-point quantum fluctuations of the 
pion field. 



This picture is similar to one which was discussed in where the relic gravitons 
production from zero-point quantum fluctuations during the cosmological expansion was 
considered. For graviton mode with momentum n the equation y" + [n^ — {R" /R)]y = 
was obtained, wheTeR{r]) is the scale factor of the metric ds"^ = R^{ri){d7f — dx^) and a 
prime represents d/dr]. One can see that this equation is equivalent to pion equation (H) 
if the scale factor R is connected with < 0^ > (t) as A(f^— < 0^ > {t)) = R"/R. 

Our problem now is to present the quantum mechanical formulation in terms of pion 
creation and annihilation operators and to get the wave function of the DCC . In the 
mean field approximation the wave function |\E' >= H - ^ IV' >j ^ is the product of the wave 
functions for each mode with momentum k and isotopical index i. Later we shall omit i. 
The equation of motion (|^) for each mode 7r(fc,t) can be derived from the lagrangian 

Lk = ^rt\k,t) -^n\k,t)n\k,t) 

n\k,t) = P + X{<(j)^>{t)-v^) (6) 
The wave function \ip >^ obeys Schrodinger equation 

4|^>a; = H,m>^, = [l'Pl + ^Q\k,t)n\k)]\^>,i (7) 
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where 7r(/c) and V-,^ = —id/dTi{k) are the quantum- mechanical coordinate and momentum 
for the mode with the spatial momentum k. One can rewrite the Hamiltonian in (|^ 
in terms of creation and annihilation operators which make it diagonal at any given 
moment. It is evident that we are interested to get the wave function in terms of creation 
and annihilation operators of ordinary pions, so we must diagonalise the Hamiltonian at 
t ^ oo, when the oscillation of the < cf? > (t) will be damped. Thus we define 

a{k) = -P^+^-CkMk)^ ^ V.-i.ikMk) 



2uj{k) y^uik) 



where uj{k) = ^l{k, oo) = yk"^ + m^. It is easy to see that the Hamiltonian is 

H, = Uk)[l + ^^]a^Ck)a{k) + + ^t2(^)] 

^ uj^{k) Auj{k) 

The a^(fc) and a^'^(k) terms in Hamiltonian is the reason that H transforms the initial 
vacuum state |0 > into a squeezed state S{r,(j))\0 > @. 

To calculate the squeezing and phase parameters r and (f) let us consider the solution 
of the Schrodinger equation (0) in the coordinate (here it is vr) representation ( we shall 
omit the label k for moment) 

^(vr, t) = C{t) exp{-B{t)Tr^) (10) 

and for B = uj/2 this wave function describes the vacuum state. For all other values 
this wave function describes the squeezed state @) where the parameters r and are 
connected with B by the relation (see also |T9| and ||21|| ) 

uj coshr + exp(2z0) sinhr 



B 



2 coshr — exp(22(/)) sinhr 



w2 + 4|5P 1 ImB , , 

cosh2r = —; — „ „ ; sin 20 = — — — - (11) 
AujReB smh 2r ReB ^ ' 



Substituting ([10| ) into (0) one gets the equation for B{f) 

W = 2B^ - ^n\t) (12) 
which means that B{t) is related to the solution of the classical equation (|])[B 

B{t) = -'-^y m+^'m{t) = o (13) 
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The last equation can be viewed as a Schrodinger equation describing the wave function 
ip{t) of a "particle" with mass m = 1/2 on a line t having energy k"^ and moving through 
the potential barrier V{t) = — A(< 0^ > (t) — v 



2^ 



- ^ + x{v'- <<P'> {t))m = em (m) 

Far from the barrier, i.e. at t ^ ±oo one has V{±oo) = —m^ and general solution of the 
Schrodinger equation aX t ^ ±cxd is the superposition of the left and right moving waves 

ijjit) = Ae-^^(*^)* + B*e+'^^^'^' t +00 

m = Ce-*^(^')* + D*e+''^('=)* t -oo (15) 



where uj{k) = \fW^^^^m?' . Due to the unitarity the total fluxes at t ^ ±oo must be equal 
\A\^ — |_Bp = |Cp — and one can find 

A = coshrC - e^^^sinhrD* 
B* = -e-^''' sinh rC + cosh rD* (16) 

where 6 is the scattering phase and factor r is defined by the probability of the transition 
through the barrier. 

Let us remember that we are starting from the vacuum at t — oo, i.e. from B = uj/2, 
so one must have C = 0. This means that at the left (large negative t) we have only left 
moving outcoming wave /)*e+*'^('^)*. At the right (large positive t) one has both left 
and right moving waves, i.e. the incoming _B*e+*'^('^')* and reflected Ae~^^^^'^* waves. The 
tunneling coefficient can be obtained from ([TBI) by puting C = 



\D 



2 



^ I -Bp cosher ^ ^ 



Now let us calculate B{t) = —{i/2){ip/ip) at large positive t. Using ([T5| ) one can find 
after simple calculations: 

^ coshr + exp[2i(^ — co'(A;)t)] sinhr 
2 coshr — exp[2i(^ — uj(k)t)] sinhr 

in complete agreement with ([Tl|), where now the phase factor cj) = 9 — uj{k)t. 

The squeezing parameter r depends on the absolute value k of the mode momentum 
k of course and as we had found it is determined by the probability of tunneling through 
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the potential barrier V(t) = A(< 0^ > (t) — t>^) . Let us note that tunnehng takes 
place precisely when — V{t) < 0, i.e. when the classical long wavelength modes are 
exponentially amplified, and we see once more that squeezing is ultimately connected with 
the exponential growth of the classical long wavelength modes. Using the quasiclassical 
approximation it is easy to calculate the squeezing parameter r{k): 



1 



cosh^ r(k) 

Thus for small k and large r{k) one has 



exp{-2Re J dt^V{t) - P) (19) 



r{k) = 2Re / dt^X{v^- < 02 > (t)) - P (20) 



and the squuezing for each mode k is determined by the function < cfp' > (t) - this is the 
only input information we must know to calculate the DCC wave function. The average 
number of particles in each mode is 

< Nk >= sinh^ rfe = exp(2i?e J dt^V{t) - F) - 1 (21) 

and we see that < Nj. > sharply decrease with the increase of k. Let us consider the 
simple model for < 0^ > (t) assuming that < 0^ > (t) = for t G (0,m~^) and then 
becomes equal to its usual v.e.v. One can estimate 



r(k) = ^fl!^, 



2^ (22) 



Of course this is the very crude estimate, however it demonstrates the qualitative features 
of the phenomenon - the sharp exponential dependence on k and large amplification dactor 
in r which is of order of mg-r, where r is the characteristic time of damping of < 0^ > (t) 
oscillations. 

Another interesting feature is the phenomenon of bunching and super-Poissonian 
statistics [jl8|, ||2^. One may consider the second-order correlation function 



2, , < iV(t)iV(0) > , , 

9\t) = ^I^^J (23) 

which gives us the relative probability to measure two particles in an interval t and g'^{0) 
measures the probability of simultaneous detection. For coherent state one has g'^{t) = 1 
which means that for coherent state the detection events are not statistically dependent. 
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For squeezed state with large r one gets|] fi'^(O) = 2+coth^ r > 3 and (7^(0) is increasing 
with increasing of k. The experimental measurement of g^{0) for pions with different k 
will be extremely interesting. 

Let us note that the picture we have considered here may have applications not only 
to DCC but to the gluon condensate too. In a recent paper it was shown that 
a high frequency standing wave in SU(2) gauge theory is unstable against decay into 
long wavelength modes which gives the mechanism for energy transfer from initial high 
momentum modes to final states with low momentum excitations. The obtained picture 
is similar to the amplification of the pion soft modes and it is possible to obtain ther 
same squeezed description of the final state for the soft gluon modes. It is interesting 
that because the squeezed final state is quadratic in field one can get gluon condensate 
<Gl^>. 

In conclusion we would like to stress that in the mean field (Hartree-Fock) approxi- 
mation the wave function of DCC is the product of the the squeezed state wave functions 
of all three pions 

1^ >Dcc= I[[]lS{r{k),m)\0 >,] (24) 
«=i k 

— * — # 

where the universal functions r{k) and (j){k) are completely determined (in the mean 
field approximation) by the only function < 0^ > (t). This wave function leads to some 
interesting predictions about k"^ dependence of the observables which will be interesting to 
check in experiment. It is also interesting to go beyond the Hartree-Fock approximation 
and get the effects of correlations between different pion modes. 

I am grateful to A. A. Anselm, D. Gross, K. Rajagopal, S.G. Matinyan B. Miiller, A.M. 
Polyakov, W.D. Walker and F. Wilczek for interesting and stimulating discussions. This 
work was supported by the National Science Foundation grant NSF PHY90-21984. 
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